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5-H ! 1. Introduction 



We use surface adapted ellipsoidal coordinates to obtain a first order post-Newtonian 
(PN) approximation to Dedekind ellipsoids with the intention of proceeding to higher 
orders. 
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Two effects may prevent the existence of non-axisymmetric, but nonetheless station- 
ary configurations in General Relativity. In the case of imperfect fluids, dissipation 
^ ' will lead to a heating of the configuration and thus to time-dependent behavior. 

This was rigorously proved in Ref. |2J However, dissipation is completely neglected 
I in perfect fluids, but the configuration can still lose energy due to gravitational ra- 

i diation. The energy loss is known from the linearized theory to be associated with 

the time variation of the Newtonian quadrupole moment. By basing one's consid- 
' erations on a Newtonian mass distribution which is non-axisymmetric, but has a 

. constant quadrupole moment, one might hope to find a fully relativistic solution, 

which is also non-axisymmetric but does not contain gravitational radiation. The 
possibility or impossibility of such configurations has not been proved yet. 

Since solving the complete problem is not feasible, a PN approximation of an 
^ ' appropriate Newtonian solution, the Dedekind ellipsoids, is performed. These ellip- 

soids are triaxial with the axes ai > 02 > 03; the matter content is described by 
a perfect fluid, whose elements move on ellipses perpendicular to the shortest axis. 
The sequence of Dedekind ellipsoids forms a one parameter family, where for ex- 
ample ^ G (0, 1) serves as the parameter. Chandrasekhar and Elbert investigated 
the deformed figures of the Dedekind ellipsoids in the PN approximation to first or- 
der.l^ These calculations were carried out in Cartesian coordinates relying heavily 
on Newtonian results obtained in Ref. [2] This approach would become increasingly 
cumbersome for higher order corrections. On the other hand, the PN approxima- 
tion of Maclaurin spheroid^?' shows that spheroidal coordinates can be implemented 
advantageously. Therefore, we attempt to do the same here and employ ellipsoidal 
coordinates, at least in a surface adapted form. 

The connection to known solutions can be made by requiring that in the limit 
of axial symmetry, ^ — !■ I, the PN approximation of the Dedekind ellipsoids be 
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equal to the PN approximation of the Maclaurin spheroids at the bifurcation point 
(1.716... as — ai — 02). This limit is thoroughly investigated in Ref. ^ for the post- 
Newtonian Dedekind ellipsoids obtained in Ref. |31 The other end of the Dedekind 
sequence (— — !■ 0) describes a static line of mass; therefore, it generates a Weyl 
field. 

2. The formalism and the equations 

To exploit the stationarity maximally, we use the projection formalism developed 
by Geroch in Ref. 5. If the time coordinate a;" is chosen to be Killing, the metric is 
decomposed into F = goo , hat = gab and = goa , where is the Killing vector and 
gfii, the metric. Greek indices run from to 3 and Latin from 1 to 3. Furthermore, 
we have a mass density p, which we assume to be a constant not subject to the PN 
approximation, a four- velocity and a pressure p. These functions are expanded in 

some dimensionless, relativistic parameter e ^ 1. The PN expansion of a quantity 

00 

/ is then given by / = ^ fM^^^ where denotes the Newtonian quantity. 

n— no 

We use surface adapted ellipsoidal coordinates. The common ellipsoidal coor- 
dinates {A,/z, !/}, as defined e.g. in Ref. 1, have the following coordinate surfaces: 
A — constant describes confocal ellipsoids, /i — constant confocal one-sheeted hy- 
perboloids and v = constant confocal two-sheeted hyperboloids. The surface of a 
Newtonian Dedekind ellipsoid is given by A = oi. Assuming we were to have already 
solved the entire problem, we can write the surface of the resulting configuration as 

A(/i, J^) = ai(l + S*^^"' (/i, i/)e^"). By introducing a new coordinate 

n=2 

00 

A' = A(l + ^^(2„)(^^^)^2„)-i^ (1) 

the surface of the configuration is still defined by X' — ai. 

In these coordinates, h^fj is the Euclidean metric in ellipsoidal coordinates and 

(2) 

the next non-vanishing order h^f^ is determined by a linear system of second order 

partial differential equations, whose solution is given by -2F(2)/i^"j^\ where F^^) 

(3) 

denotes the Newtonian potential. The equations for can in an appropriate gauge 
be reduced to a Poisson equation where the inhomogeneity is given by the Newtonian 
velocity distribution. This equation can be solved by separation of variables and an 
expansion in ellipsoidal surface harmonics.'^ Therefore, and h^^jj are obtained 
without using the PN corrections of the surface, the four-velocity and the pressure. 
The structure of the equation for F^^^ is given by 

where the function gi depends on Newtonian quantities Z^""-*, on the coordinates 
and linearly on S'^^^; A denotes the Laplace operator in Euclidean space. If 5(2) IS 
expanded in ellipsoidal surface harmonics, Eq. ([2]) can be solved and the solution is 
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determined by the expansion coefficients of 5^^^. The boundary conditions for Eq. 
([2]) are given by the assumption of asymptotic flatness. The functions describing 
the matter and F^^^ can be further restricted using the Bianchi identity 



where gs is again hnear in the PN quantities. These equations have to be consid- 
ered in the matter configuration, i.e. in the Newtonian efiipsoids, because of our 
surface adapted coordinates. The integrability condition of Eq. (|3b|l yields a equa- 
tion relating S*'^) and u^^^ provided that the general solution for F^"') is inserted. 
Furthermore, Eq. (j3b|) determines the PN correction of the pressure up to a constant. 
Eq. (jSaJ restricts the four-velocity, but does not determine it completely. There is 
still the freedom to add a source-free vector field to the spatial components of the 
four-velocity. There are two more conditions, which are satisfied by a stationary 
solution. The pressure vanishes at the surface as well as the component of the four- 
velocity normal to the surface. These conditions take on a particularly simple form 
in the coordinate system mentioned above since they reduce to p{X' = ai, ^,,v) = 
and (A' — ai,^,^) — 0. Solving these equations to the first PN order, the four- 
velocity and the surface corrections are further fixed. 
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